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An XMM-Newton observation of Ton S180: Constraints on the 
continuum emission in ultrasoft Seyfert galaxies 
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ABSTRACT 

We present an XMM-Newton observation of the bright, narrow-line, ultrasoft Seyfert 1 galaxy 
Ton S180. The 0.3-10 keV X-ray spectrum is steep and curved, showing a steep slope above 
2.5 keV (F ~ 2.3) and a smooth, featureless excess of emission at lower energies. The spec- 
trum can be adequately parameterised using a simple double power-law model. The source is 
strongly variable over the course of the observation but shows only weak spectral variability, 
with the fractional variability amplitude remaining approximately constant over more than a 
decade in energy. The curved continuum shape and weak spectral variability are discussed 
in terms of various physical models for the soft X-ray excess emission, including reflection 
off the surface of an ionised accretion disc, inverse-Compton scattering of soft disc photons 
by thermal electrons, and Comptonisation by electrons with a hybrid thermal/non-thermal 
distribution. We emphasise the possibility that the strong soft excess may be produced by 
dissipation of accretion energy in the hot, upper atmosphere of the putative accretion disc. 

Key words: galaxies: active - galaxies: Seyfert: general - galaxies: individual: Ton S180 - 
X-ray: galaxies 



1 INTRODUCTION 

Many Seyfert 1 galaxies and radio-quiet quasars possess X-ray 
spectra that steepen below ~ 2 keV. This 'soft excess' emission, 
usually defined as the excess emission over and above the hard 
X-ray power-law, was first seen in spectra obtained by HEAO-1 
A2 (Pravdo et al. 1981), EXOSAT (Arnaud et al. 1985; Turner & 
Pounds 1989) and Einstein (Bechtold et al. 1987). Its physical ori- 
gin remains uncertain, as does its connection with the so-called 'big 
blue bump' emission rising through the ultraviolet, although it is 
often associated with thermal emission from the putative accretion 
disc (e.g. Arnaud et al. 1985; Czerny & Elvis 1987). 

Boiler, Brandt & Fink (1996) and Laor et al. (1997) showed 
that the objects with the steepest soft X-ray spectra (implying 
the strongest soft excesses) tend to have relatively narrow optical 
H/3 lines; many are classified optically as narrow-line Seyfert Is 
(NLSls; Osterbrock & Pogge 1985). These 'ultrasoft' Seyferts of- 
ten show other notable properties such as very rapid X-ray variabil- 
ity (e.g. Forster & Halpem 1996; Boiler et al. 1997; Turner et al. 
1999; Leighly 1999a; Brandt et al. 1999) and strong optical Fe II 
emission (e.g. Lawrence et al. 1997; Vaughan et al. 2001). The X- 
ray spectral form. X-ray variability and optical broad-line proper- 
ties seem interconnected. 

Here we report the results of a 30 ksec XMM-Newton observa- 
tion of Tonantzintla S 1 80 (2: — 0.062) obtained as part of a guaran- 
teed time programme (PI: Th. Boiler) to study the timing and spec- 



tral properties of ultrasoft NLSls using XMM-Newton. Ton SI 80 is 
one of the X-ray brightest ultrasoft Seyferts (Vaughan et al. 1999; 
Leighly 1999a) and its luminosity is such that it is often classi- 
fied as a quasar. Previous X-ray observations with ASCA (Turner, 
George & Nandra 1998; Vaughan et al. 1999; Leighly 1999b; Bal- 
lantyne, Iwasawa & Fabian 2001) and BeppoSAX (Comastri et al. 
1998) showed a strong excess of soft emission and tentative ev- 
idence for an ionised iron line. The more recent high-resolution 
Chandra LETGS spectrum (Turner et al. 2001b) showed no strong, 
narrow absorption or emission features and showed the soft X-ray 
continuum is smooth and featureless. 

The rest of the paper is organised as follows. In the follow- 
ing section the observation and basic data reduction procedures are 
outlined. Section ^ describes the X-ray spectral fitting results. This 
is followed by an analysis of the variability properties of Ton S 1 80 
in section Q The implications of these results are discussed in sec- 
tion |[ 



2 OBSERVATION AND DATA REDUCTION 

XMM-Newton (Jansen et al. 2001) carries three European Pho- 
ton Imaging Cameras (EPIC), specifically two MOS (Turner et al. 
2001a) and one pn (Struder et al. 2001) CCD camera. In addition 
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there are two Reflection Grating Spectrometers (RGS; den Herder 
et al. 2001) and an Optical Monitor (OM; Mason et al. 2001). 

XMM-Newton observed Ton SI 80 on 2000 December 14 (rev- 
olution 0186) for a duration of ~30 ksec, during which all in- 
struments were operating nominally. The MOSl camera was in 
timing mode, which severely complicates the analysis, and these 
data are ignored in the present paper. Both the M0S2 and pn 
were operated in small window mode to reduce "pile up" from 
this relatively bright (few ct s^^) source, and all EPIC cameras 
used the medium filter. The RGS was operated in standard (Spec- 
tro+Q) mode. Extraction of science products from the Observation 
Data Files (ODFs) followed standard procedures using the XMM- 
Newton Science Analysis System v5.3 (SAS). 

The raw data were processed to produce calibrated event lists 
and screened to remove unwanted hot, dead or flickering pixels, 
and events due to electronic noise. Light curves extracted from 
these event lists showed the background to be stable throughout 
the duration of the observation. The total amount of "good" expo- 
sure time selected was 28, 644 s and 20, 338 s for M0S2 and pn, 
respectively, and 30, 232 s and 29, 321 s for RGS 1 and RGS2. (The 
lower pn exposure is due to the lower "live time" of the pn camera 
in small-window mode, ~71 per cent; Striider et al. 2001). 

Source data were extracted from circular regions of radius 
35" from the MOS2 and pn, and show no signs of pile-up. Back- 
ground events were extracted from off-source regions. Events cor- 
responding to patterns 0-12 (single-quadruple pixel events) were 
extracted from M0S2 and patterns 0-4 (singles and doubles) were 
used for the pn analysis, after checking for consistency with the 
data extracted using only single pixel events (pattern 0). Stan- 
dard redistribution matrices (m2_r6_all_15 . rmf for MOS2 and 
epn_sw20_sdY9 . rmf for the pn) were used, and ancillary re- 
sponse files were generated using ARFGEN V 1.48. 8. First-order 
RGS spectra were extracted using the RGSPROC V 1 . 3 . 3 script, and 
appropriate response matrices were generated using RGSRMFGEN 
V 1.44. 5. The total number of "good" source events extracted was 
2.9 X 10^ for the pn, 8.8 x 10* for MOS2 and 1.4 x lO" for each 
RGS. 



3 SPECTRAL ANALYSIS 

The source spectra were grouped such that each spectral bin con- 
tains at least 20 counts, and they were fitted using the XSPEC V 1 1 . 1 
software package (Arnaud 1996). The quoted errors on the derived 
best-fitting model parameters correspond to a 90 per cent confi- 
dence level for one interesting parameter (i.e., a Ax^ = 2.7 crite- 
rion) unless otherwise stated. Values of Ha = 70 km s^^ Mpc~^ 
and go = 0.5 are assumed throughout the paper, and fit parameters 
are quoted for the rest frame of the source. 

3.1 The 2.5-10 keV spectral form 

The two EPIC spectra were first examined separately to test for 
discrepancies. The data were fitted over 2.5-10 keV using a model 
comprised of a power-law modified by Galactic absorption (with 
fixed column density Nh = 1-55 x 10^" cm~^; Dickey & 
Lockman 1990). This provided a reasonable fit to both spectra 
(X^ = 115.0/153 degrees of freedom and = 339.8/336 dof 
for the MOS2 and pn data, respectively) with consistent normal- 
isations and photon indices (T = 2.29 ± 0.10 for MOS2 and 
r = 2.25 ± 0.06 for pn). 

Given this good agreement between the two spectra over the 
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Figure 1. EPIC pn (upper data) and MOS2 spectra in the 2.5-10 keV range 
compared to a simple power-law model (the energies have been shifted into 
the source frame). The upper panel shows the count spectra (crosses) and 
the power-law model folded through the detector responses (histograms). 
The lower panel shows the fit residuals. 

2.5-10 keV energy range, the data were fitted simultaneously in 
the rest of the spectral analysis, but in each fit the residuals from 
the two spectra were examined separately to check for inconsisten- 
cies. The relative normalisations were left free to allow for slight 
differences in absolute flux calibration between the two detectors. 
Figure |^ shows the two EPIC spectra compared to a simple power- 
law fitted over 2.5-10 keV. As expected from the fits to the individ- 
ual spectra, this model (F — 2.26lg ^2) excellent agreement 
with the data (x^ = 455.1/490 dof). 



3.1.1 Iron Ka emission 

The pn spectrum shows positive residuals around ~ 7 keV which 
could indicate spectral features from iron K. To test for this an in- 
trinsically narrow {g = 1 eV) Gaussian emission line was added 
to the model. With the narrow line energy fixed at 6.4 keV the fit 
was not significantly improved; the 90 per cent upper limit of the 
equivalent width of a narrow, neutral iron line is EW < 60 eV. 
Allowing the line energy to be free provided an improvement to 
the fit (x^ = 447.3/488 dof) but the best-fitting line energy of 
E = 7.01 ± 0.06 keV is marginally too high for Ka emission 
from even Hydrogen-like iron. The equivalent width of this line 
is EW = 72 ± 43 eV. According to an F-test, upon adding 
two parameters (line flux and energy) the fit is significantly im- 
proved at 98.5 per cent confidence.. Allowing the line width to 
be free improved the fit further (x^ ~ 441.5/487 dof) with 
a best-fitting width of ct = 0.53^0 33 keV and an energy of 
E — 7.01 ± 0.31 keV, which is consistent with Ka emission from 
H-like iron. An i^-test suggests the fit is improved, at 98.8 per cent 
confidence, upon allowing the line width to be free. This broad line 
is strong, with an equivalent width of EW = 282t^g3 eV. The 
high energy derived from the line might instead be indicating the 
presence of iron K absorption (at i5 > 7.1 keV), but this seems 
unlikely. Replacing the broad Gaussian with an absorption edge 
gave a worse fit (x^ = 451.8/488 dof), with edge parameters 
E = 8.33t[; j^ keV and r = 0.15 ± 0.14. 

Using a DISKLINE line profile (Fabian et al. 1989) to model 
the residuals around 7 keV provided a comparable fit to the broad 
Gaussian (x^ ~ 444.2/487). An F-test suggests this improve- 
ment to the fit, compared to the simple power-law model, is sig- 
nificant at 99 per cent confidence. The best-fitting parameters of 
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the line were the following: rest-frame energy E = 6.48lQ Qg keV, 
EW = 5231:111 eV and inclination angle i = 65tf deg. Allow- 
ing the other free parameters (inner radius ri„, outer radius rout, 
and disc emissivity index q) to be free did not improve the fit sig- 
nificantly, and thus these parameters were kept fixed at rin = 6rg, 
Tout = lOOOrg and g = 3. As can be seen from figure |l|, the ex- 
cess at ~ 7 keV in the pn data is not obvious in the M0S2 data; 
however, the two spectra are formally consistent with one another. 
When the normalisation of the DIS KLINE was allowed to be free be- 
tween the two detectors, the best fitting normalisations were consis- 
tent, although the normalisation form the M0S2 spectrum is poorly 
constrained and is consistent with zero. 

The above results are suggestive of "reflection" off a (pos- 
sibly ionised) accretion disc. To provide a more physically real- 
istic description of emission from an ionised disc, the model of 
Ross & Fabian (1993) (see also Ballantyne et al. 2001) was fitted 
to the data (assuming solar abundances). The disc reflection spec- 
trum was blurred using the LAOR kernel (Laor 1991) to simulate 
Doppler and gravitational effects around a black hole. The best fit- 
ting ionised disc model gives a good fit (x^ = 445.0/487 dof) 
with r = 2.43 ± 0.13, R = 1.5^1:1 and fairly low ionisation 
(log(^) < 1.9). The LAOR profile parameters were only poorly 
constrained and were kept fixed (as above). 

We conclude that the presence of a broad, ionised iron line 
(as previously suggested by Comastri et al. 1998 and Turner et al. 
1998) is consistent with the XMM-Newton data, but the line param- 
eters are only poorly constrained. 



3.2 The broad-band X-ray spectrum 

Having parameterised the spectrum above 2.5 keV the analysis was 
extended to cover the full useful spectral bandpass of XMM-Newton 
by including the EPIC data down to 0.3 keV and the RGS data. 
However, there are significant differences between the MOS2 and 
pn spectra below 1 keV; the difference is most pronounced around 
0.56 keV. There is a similar feature present in the calibration obser- 
vation of Mrk 421, which is almost certainly related to the instru- 
mental O K-edge at 0.537 keV and seems to be most pronounced 
in data taken in pn small-window mode (as here). The residuals in 
the MOS camera are at a much lower level (<J 5 per cent) in this 
spectral region. For these reasons a conservative restriction was ap- 
plied to the pn data, the spectrum was fitted between 1.2-10 keV 
only, while the M0S2 data were used over the 0.3-10 keV band. 

An extrapolation of the simple 2.5-10 keV power-law model 
down to 0.3 keV (Figure ^ shows a smooth upturn in the EPIC 
spectrum: the soft excess. There does not seem to be an obvious 
transition from hard power-law to soft excess in the spectrum, on 
the contrary, the spectrum is gently curving up to at least 2 keV. 

The average observed flux in the 0.3-10 keV band was 2.2 x 
10~^^ erg s~^ cm~^ and the corresponding (unabsorbed) luminos- 
ity was ~ 2.5 X lO** erg s~^. The bolometric luminosity is most 
likely > 10'*^ erg s~^ (this is a conservative estimate as Ton SI 80 
probably harbours a strong EUV excess; Turner et al. 2002). As- 
suming Ton SI 80 is radiating at the Eddington limit (or lower) then 
it requires a black hole mass Al > 10^ Mq (see also Turner et al. 
2001b). 



3.2.1 RGS spectrum 

The 0.35-1.5 keV RGS data confirm that the soft excess emission 
in Ton SI 80 is a smooth continuum, as seen by Turner et al. (2001b) 
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Figure 2. EPIC pn (upper data) and MOS2 spectra. The upper panel shows 
the count spectra, and the lower panels show the ratio of data to three dif- 
ferent models. The second panel shows the results of fitting a single power- 
law model only to the data above 2.5 keV (as in figure |l|) and extrapolating 
to lower energies, revealing the curved continuum. The bottom two panels 
show the results of fitting with a power-law plus blackbody or a two power- 
law model. (The data were rebinned for display purposes.) 



from Chandra LETGS data. Figure |3| shows the background sub- 
tracted RGS spectrum of Ton SI 80 in flux units (RGSl and RGS2 
combined). An absorbed power-law provides a good fit to the data 



1315.2/1256 do/) with the absorbing column fixed at the 



Galactic value and a power-law slope of F = 2.77 ± 0.03. There 
are no significant sharp residuals in the RGS spectrum, the only ob- 
vious deviation from the power-law model is around 0.5-0.6 keV 
and is probably associated with the O K-edge in the detector re- 
sponse. The equivalent widths of any narrow absorption features 
are restricted to EW < 1.0 eV (90 per cent upper limit) in the 
0.35-0.70 keV energy range, EW < 1.7 eV in the 0.7-1.0 keV 
range and EW < 10 eV in the 1.0-1.5 keV range. Thus the soft 
excess in Ton S 180 is a smooth continuum resembling a power-law. 
In the following sections the EPIC spectra are used to test various 
physical models for the soft excess emission. 

3.2.2 Phenomenological model 

Given that the RGS spectrum can be fitted with a steep power law, 
the broad-band EPIC spectrum was first compared to a simple dou- 
ble power-law model, similar to that used by Comastri et al. (1998) 
in their analysis of a BeppoSAX observation of Ton SI 80. Fitting 
with a broken power-law model gave = 919.1/898 dof with 
slopes Fi = 2.98 ± 0.02, r2 = 2.22 ± 0.10 and a break en- 
ergy Sbr = 2.29 ± 0.04 keV. A slightly better fit was obtained 
using two separate power-laws (x^ — 910.9/898 dof) with slopes 
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Figure 3. Combined first order RGS spectrum (in "fluxed" units) and best- 
fitting absorbed power-law model (solid line). The lower panel shows the fit 
residuals. 
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Figure 4. Spectral energy distribution of Ton S 1 80, corrected for Galactic 
absorption. The solid curve shows the XMM-Newton spectrum (using the 
broken power-law parameteiisation discussed in section 3.2.2). The dots 



show the multiwavelength data from Turner et al. (2002) corrected to match 
the slightly different cosmology used here. 



Fi = 3.11 ± 0.04 and r2 = 1.49 ± 0.15 (see figure ^. In this 
model the softer of the two power-laws dominates the spectrum up 
to ~ 4 keV. The harder power-law is considerably flatter than in 
the broken power-law parameterisation, due to the addition of the 
steeper power-law over the whole band, and is perhaps unrealisti- 
cally flat. These simple fits to the EPIC data, and the RGS results, 
show that the soft excess has the form of a steep power-law. 

Figure places the XMM-Newton spectrum in the context of 
the multiwavelength spectral energy distribution (SED) compiled 
by Turner et al. (2002). The multiwavelength data were taken ap- 
proximately one year before the XMM-Newton observation, but it 
is clear that the X-ray spectral form did not change drastically be- 
tween observations. 



3.2.3 Ionised reflection 

The possible presence of iron Ka emission (section |3.1.l| ) may in- 
dicate significant reflection off the surface of the accretion disc. 
The reflection spectrum can be strong below ~ 1 keV and could 
contribute significantly to the soft excess. To test whether the soft 
excess in Ton SI 80 could be produced by reprocessing of pri mary 



X-rays in the disc the reflection model discussed in section 3.1.1 
was fitted to the 0.3-10 keV spectrum. The best fitting model ob- 
tained for the data above 2.5 keV gave a very poor fit to the broad- 
band data (x^ ~ 8351.1/897 dof). The ionisation state of the 
reflector means it produces strong O VII recombination emission, 
not obvious in the data, while the total reflected flux below 2 keV 
falls far short of the observed soft excess. 

Refitting this model to the full-band EPIC data still provided 
a rather poor fit to the data (x^ = 1030.4/896). The best fitting 
parameters were T = 2.58 ± 0.03, R = 1.121!^;?^, Iog(5) = 
3.71 ± 0.06 and parameters of the LAOR kernel were nn < 2.5rg 
and i = 46 ± 8° (rout was kept fixed at 400r-g). In this model 
the disc surface is highly ionised and produces a smooth excess of 
emission in the soft band, partly from bremsstrahlung in the irradi- 
ated skin of the disc, with relatively weak spectral features. How- 
ever, this model leaves significant residuals in the fit above 7 keV, 
where it predicts a strong, ionised iron edge. Thus the low-energy 
data seemed to require a highly ionised reflector, while the higher- 
energy spectrum suggested a reflector with much lower ionisation. 

Allowing for two reflection components, with different ion- 
isation parameters, improved the fit (x^ = 964.7/894 dof). In 
this fit the continuum slope was F = 2.61 ± 0.02 incident on 
reflectors with ionisation parameters log(^i) = 4.4 ± 0.3 and 
log(^2) = 1.8±0.1, the LAOR profile parameters were Tin < l-9rg 
and i — 41 ± 6°. In this model the highly ionised reflector had 
relative strength Ri « 0.5 while the colder reflector was much 
stronger, with R2 ~ 2.4. Both reflection components contributed 
to the soft excess emission, but the strong colder reflector improved 
the fit to the hard X-ray data, although significant residuals still re- 
mained above 7 keV. The small inner radius of the LAOR kernel was 
required in order to smooth sufficiently the recombination features 
present in the reflection spectra. 



3.2.4 Purely thermal soft excess 

A power-law plus single blackbody component, to model the soft 
excesSjprovided a poor fit to the data (x^ ~ 1123.6/898 dof, see 
figure |2). Including a second blackbody component significantly 
improved the fit (x^ = 912.2/896 dof), while a power-law plus 
three blackbodies provided an excellent fit (x^ ~ 887.8/894 dof) 
with the following parameters: F = 2.24 ± 0.02, kTi = 33tf eV, 
kT2 = 103 ± 5 eV, fcTa = 248 ± 15 eV. This multiple blackbody 
parameterisation is similar to that used to characterise the broad- 
band X-ray spectra of other ultrasoft Seyferts, e.g., RE J1034-I-396 
(Pounds, Done & Osborne, 1995), PKS 0558-504 (O'Brien et al. 
2001) and IH 0419-577 (Page et al. 2002). 

However, while this multiple blackbody model did provide an 
adequate fit to the data, the temperatures and luminosities of these 
components are at odds with those expected for a simple accretion 
disc. A standard accretion disc around a lO'^ Mq black hole should 
have a peak temperature <. 40 eV, much lower than the hottest 
temperature in the multiple blackbody model. The observed tem- 
peratures and luminosities predict sizes of the blackbody emitting 
regions of rae ~ 3 x 10^" cm for the hottest and t-bb ~ 10^"^ cm 
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for the coolest, which are also difficult to reconcile with a simple 
accretion disc around a lO'^ M© black hole (r-g ~ 1.5 x 10^^ cm). 

Fitting the data with a more realistic thermal accretion disc 
spectrum, using the DISKPN code (Gierlihski et al. 1999), provided 
a poor fit to the data (x^ = 1073.7/898 dof), because the model 
spectrum is not broad enough, and again it gave an unreasonably 
high temperature (kT — 131 ± 5 eV). The strong, rapid variabil- 
ity of the soft excess (section ^ also seems incompatible with the 
expected timescales in a simple accretion disc (see also Boiler et 
al. 1997). These results would seem to reject an origin for the soft 
excess in terms of (un-modified) thermal emission from an accre- 
tion disc. But, as discussed by e.g. Bechtold et al. (1987), Czerny 
& Elvis (1987) and Ross, Fabian & Mineshige (1992), inverse- 
Compton scattering in the upper layers of the disc atmosphere 
should significantly affect the observed high-energy tail of the disc 
emission. This possibility is discussed below. 



3.2.5 Thermal Comptonisation 

Inverse-Compton scattering of soft photons in an optically thick 
medium (i.e., r > 1; where the diffusion approximation is valid) 
produces an approximately power-law spectrum at low energies. 
Thus, as the soft excess in Ton SI 80 resembles a steep power- 
law, a Comptonised blackbody component should provide a rea- 
sonable fit. The energy index of the Comptonised spectrum is 
determined by the Compton y-parameter (j/ ~ 4Gr^, where 
© = kT /meC? is the dimensionless electron temperature) ac- 
cording to Q « y'9/4 + 4/y — 3/2 (e.g., section 7.7 of Rybicki 
& Lightman, 1979). This means that it is not possible to solve 
for both r and kT from knowledge of the spectral slope alone. 
Usin g the double power-law parameterisation of Ton SI 80 (sec- 
tion 3.2.21 gave a slope of a = 2.1 for the soft excess, which 
corresponds to j/ ~ 0.37 and it is possible to obtain this value. 



and hence the observed spectral slope, for a range of r and kT 
values. A similar situation occurs when considering the optically 
thin regime, in which case the spectral slope is determined by 
a - In r/ In (1 49 + 166^). See Zdziarski (1985) and refer- 
ences therein for more details. 

The CompTT code (Titarchuk 1994) was used to model 
Comptonisation of soft photons in a thermal plasma. A power- 
law plus Comptonisation component gave a good fit to the data 
(X^ = 903.7/896 dof) with V = 1.64 ± 0.18 and the soft 
excess modelled by Comptonisation of a thermal spectrum from 
a fcTbb = 60 ± 7 eV source. However, for the reasons men- 
tioned above, the temperature and optical depth of the Compton- 
ising plasma are strongly covariant parameters, and thus cannot be 
constrained simultaneously. 

Fitting the spectrum with just two Comptonisation compo- 
nents, similar to the model used by O'Brien et al. (2001) and Page 
et al. (2002), gave a good fit (x^ = 891.9/895 dof). The seed pho- 
ton temperature was fcTbb = sV, but the parameters of the 
two plasmas were poorly constrained due to the degeneracy men- 
tioned above. (Similar fits were obtained using the EQPAIR code, 
discussed below, to model two purely thermal plasmas.) It seems 
reasonable to associate the source of the soft seed photons with the 
thermal emission from the inner regions of an optically thick ac- 
cretion disc. The existence of two discrete Comptonising plasmas 
is not so easy to explain, especially given that they varied almost 
identically on short timescales (section M). 



3.2.6 Non-thermal Comptonisation 

The double Comptonisation model described above has the two 
spectrally identified continua (soft and hard power-laws) originat- 
ing in two distinct thermal Comptonising plasmas. An alternative 
is that the whole spectrum is produced by a single plasma with a 
hybrid thermal/non-thermal electron distribution (see Coppi 1999 
for a discussion of the physics of such plasmas). Gierlihski et al. 
(1999) modelled the spectrum of the Galactic Black Hole Candi- 
date (GBHC) Cygnus X-1 in its high/soft state in terms of a hy- 
brid plasma. In this model the Cyg X-1 spectrum is comprised 
of a blackbody from the accretion disc, plus a component due to 
inverse-Compton scattering of the disc photons in a corona contain- 
ing both thermal and non-thermal electrons. The non-thermal tail to 
the electron distribution produces the steep power-law tail of the X- 
ray spectrum dominating above 10 keV (and observed up to at least 
~ 800 keV) and the thermal electrons produce a "soft excess" at a 
few keV between the non-thermal tail and the accretion disc spec- 
trum. The definite detection of a non-thermal Comptonisation com- 
ponent requires high-energy data, which are not available here, but, 
as first noted by Pounds et al. (1995), the X-ray spectra of ultrasoft 
Seyferts do resemble that of Cyg X-1 in its high/soft state. Thus 
it seems reasonable to test whether a hybrid thermal/non-thermal 
plasma is a viable model for the X-ray continuum of Ton SI 80. 

The EQPAIR code (Coppi 1992; Gierlihski et al. 1999) was 
used to compute the spectrum formed by Comptonisation of a 
blackbody photons in a plasma with an hybrid electron distribu- 
tion. Non-thermal (relativistic) electrons are injected (by some un- 
specified process) as a power-law in Lorentz factor, 7, into a back- 
ground thermal plasma, the rate of injection is given by: TVinj (7) oc 
-y^Tinj (ijgt^yggu -yjj^jj^ and 7max). The total optical depth (includ- 
ing that produced by e^ pairs) and the equilibrium temperature 
of the thermal component (after accounting for Coulomb heating 
by non-thermal e* and Compton heating/cooling) are computed 
self-consistently once the initial parameters have been specified. 
For a detailed discussion of the processes accounted for by this 
code see Gierlihski et al. (1999) and Coppi (1999). The impor- 
tant parameters of the model are the following: the temperature 
of the input thermal disc spectrum fcTbb; the soft compactness 
Is — LbO-t /rnieC^ , which corresponds to the power supplied by 
soft seed photons (Ls is the total seed photon luminosity and r is 
the radius of the emission region); the similarly defined thermal 
compactness Ith, corresponding to the heating supplied to thermal 
electrons; the non- thermal compactness Inth, which gives the rate 
of injection of non-thermal electrons; the index of the non-thermal 
electron injection spectrum Finj; and the optical depth Tp of the 
background electron-proton plasma. 

With so many parameters (some of which were covariant in 
the fit) the x^ -space contained many local minima. In order to re- 
duce the complexity of the fit some parameters were fixed at phys- 
ically reasonable values (these are discussed further below). The 
details of the fitted model parameters will obviously depend on the 
chosen values for those parameters kept fixed. Thus, while the fol- 
lowing fits demonstrate that the EQPAIR model can reproduce the 
shape of the XMM-Newton spectrum, the exact values of the fitted 
model parameters should be treated with some caution. 

The total compactness I — Is + Ih (where = ^th + ^nth is 
the total hard compactness) is not constrained by the XMM-Newton 
spectra. If the spectrum extends out to 500 keV a high compact- 
ness will lead to efficient pair production, but the strongest observ- 
able consequence of the presence of pairs is an annihilation line at 
511 keV. In the absence of this information the soft compactness 
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was therefore kept fixed at Ib ~ 30 (the lower Umit obtained in sec- 
tion^) and the ratios Ih/la and ?nth/'h were used as free parameters 
in the fitting. Assuming the corona covers the inner disc, the former 
is related to the ratio of power dissipated in the corona compared to 
the disc, and the latter gives the fraction of power deposited in the 
corona that goes into the non-thermal component. The radius of the 
emitting region was kept fixed at r = 10^'* cm, this is < 70rg for 
a ^ 10^ M0 black hole, and at the above assumed value of Zs this 
gives a reasonable soft X-ray luminosity. The temperature of the 
seed blackbody was fixed at fcTbb = 50 eV and the minimum and 
maximum Lorentz factors for the non-thermal injection spectrum 
were 7miii = 1-3 and 7max = 1000 (the resulting model spectrum 
is largely insensitive to the chosen value of 7max when Finj > 2). 

The EQPAIR model gave a good fit (x^ = 892.9/897 dof) 
with the following parameters: /h/^s = 0.74 ± 0.03, Znth/^h = 



0.78 ± 0.02, 



7 0+"-* 

' -"-1.8 



and Fin 



2.84 



+0.26 



Thus there 



is marginally more power supplied by the seed photons than the 
coronal electrons, and most of the electron heating is in the form 
of non-thermal electrons which are injected with a steep spectrum. 
The soft excess is produced by Compton scattering of the accretion 
disc spectrum by thermalised electrons, and the emission above 
a few keV is produced by the non-thermal component. These re- 
sults are similar to those obtained by Gierliiiski et al. (1999) for 
the high/soft state spectrum of Cyg X- 1 ; the significant difference 
is that in this case the optical depth of the background plasma is 
much higher than in Cyg X-1. If the optical depth is kept fixed 
at Tp = la slightly poorer (but still acceptable) fit is obtained 
(X^ ~ 910.0/898 dof). Again the soft excess is produced by 
Compton scattering of seed photons by thermal electrons, and the 
harder X-ray emission is produced by the non-thermal electrons. 



4 TIMING ANALYSIS 

This section describes timing analyses of the EPIC data. Results 
are presented only for the pn data, which have the highest signal- 
to-noise, but the results were checked against the MOS2 data and 
the two were found to be entirely consistent. The (Icr) errors on the 
light curves were calculated using counting statistics. 

The topmost panel of figure ^ shows the broad-band (0.2- 
10 keV) EPIC pn light curve of Ton SI 80 in 200 s time bins. The 
light curve is perfectly evenly sampled, uninterrupted and shows 
a peak-to-peak variation of ~ 50 per cent over the course of 
the observation. The fastest continuous rise in the light curve oc- 
curred at around 2.4 x lO^s, where the flux increased by ~ 30 
per cent in 2000 s, corresponding to an increase in luminosity of 
AL/At « 3.8 X 10*'' erg s"^ in the 0.3-10 keV band. This was 
used to estimate the radiative efficiency of the source, assuming 
photon diffusion through a spherical mass of accreting matter, fol- 
lowing Fabian (1979): AL/At ^ 77 • 2.1 x lO"^ erg s"^. (But see 
also Brandt et al. 1999, for a discussion of the assumptions used in 
this calculation.) This gave a limit of 77 > 1.8 per cent. The rapid 
variability was also used to place a limit on the compactness pa- 
rameter, following Done & Fabian (1989). Assuming that r — cAt 
gives AL/At < I ■ 10^^ erg s"^. For Ton S180 this gave / > 34. 

Figure ^ also shows the pn light curves extracted from five 
different energy bands in 500 s bins. The light curves all show sim- 
ilar trends, albeit with relatively poor signal-to-noise in the harder 
bands. The hardness ratios of these light curves were examined as 
a first test of spectral variability. The ten hardness ratios (six of 
which are shown in figure ^ from the five light curves were each 
compared to a constant using a test. Four of the ten hardness ra- 
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Figure 5. EPIC pn light curves normalised to unit count rate. The topmost 
panel shows the full-band 0.2-10.0 keV light curve (in 200 s bins) and the 
lower panels show the light curves from various energy bands (in 500 s 
bins). The vertical bar indicates the (lineal') j/-axis scale. The mean count 
rates for the six light curves (from top to bottom) are as follows: 9.26, 3.26, 
3.53, 1.69, 0.56 and 0.22 ct s^^. 
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Figure 6. Hardness ratios for various EPIC pn light curves using 500 s time 
bins. 



tios examined were inconsistent with a constant hypothesis at > 99 
per cent confidence. These ratios, namely (0.4 — 0.8) / (0.2 — 0.4), 
(0.8-1.5)/(0.2-0.4), (1.5-3)/(0.2-0.4) and (1.5-3)/(0.4- 
0.8), are also among the highest signal-to-noise ratios. Using the 
Spearman rank-order correlation coefficient and the Kendall r- 
statistic (see e.g.. Press et al. 1992) to test for correlations, none of 
the hardness ratios was correlated with the broad-band count rate 
at greater than 90 per cent confidence. 

Figure ^ shows a comparison of pn spectra extracted from in- 
tervals of high and low flux. These were extracted from only those 
intervals when the source count rate was above 10 ct s~^ (high) or 
below 8.8 ct s~^ (low). The high and low flux spectra contained 
9.6 x 10'* ct and 8.0 x 10'' ct, respectively. The ratio of high to 
low flux spectra is consistent with a constant (x^ ~ 81.6/84 dof), 
again indicating a lack of flux-correlated spectral variability. 

The normalised RMS spectrum was calculated in order to 
quantify the degree of variability in the different energy bands. 
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Figure 7. Comparison of high and low flux spectra. The top panel shows 
the two count spectra taken from high flux (upper crosses, black) and low 
flux (lower crosses, red) intervals. The bottom panel shows the ratio of high 
to low flux spectra, which is consistent with a constant (dashed hne). 
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Figure 8. Normalised RMS spectrum based on EPIC pn light curves. The 
^var values were calculated using 500 s binned data. 



again using the 500 s binned light curves. The fractional variabil- 
ity amplitude was measured using the i<Var statistic (Edelson et al. 
2002) which is equal to the square root of the normalised excess 
variance CTxs (Nandra et al. 1997). The errors were estimated us- 
ing the formula in Edelson et al. (2002) and should be considered 
as conservative estimates on the true uncertainties. Figure ^ shows 
the fractional variability amplitude calculated in ten energy bands, 
which changes by only 3 per cent over the XMM-Newton band-pass. 
This remarkably flat RMS spectrum contrasts with the situation in 
many 'normal' spectrum Seyfert Is, which tend to show larger frac- 
tional variations at lower energies (e.g. Nandra et al. 1997). Alter- 
ing the sizes of the time bins or the energy ranges used in the light 
cim'es made little difference to the RMS spectrum since most of 
the variability power is on timescales > 500 s. 

Having established that the light curves show similar trends 
and almost identical amplitudes, the temporal cross-correlation 
functions (CCFs) were calculated in order to search for leads or lags 
between different energy bands. The Discrete Correlation Function 
(DCF) of Edelson & Krolik (1988) was used as an estimator of 



-10* 



-5000 







Time shift (sec) 



5000 



10* 



Figure 9. Temporal cross-correlation function for the 0.2-0.3 keV and 0.8- 
1 .0 keV EPIC pn light curves (using 500 s bins), showing a peak at zero lag. 



the CCF. The other nine energy bands used in the RMS spectrum 
were compared to the softest band (0.2-0.3 keV). Of the nine CCFs 
examined eight peaked at zero-lag (e.g., see figure ^ and the one 
exception (0.2 — 0.3 keV against 1.5 — 2.0 keV) peaked at a lag 
of only one time bin (500 s, with the harder band lagging the softer 
band) which we do not consider to be significant. In all cases the 
CCF was reasonably symmetric. These results indicate there is no 
evidence for a lag or lead between the various energy bands to a 
limit of ±500 s. 



5 DISCUSSION 

This paper presents an uninterrupted ~ 30 ksec XMM-Newton ob- 
servation of the ultrasoft, narrow-line Seyfert 1 galaxy Ton SI 80. 
The spectrum above 2.5 keV is slightly steeper than the norm for 
Seyfert 1 galaxies (F « 2.26 compared to 1.9) but the iron Ka line 
emission is only poorly constrained. At lower energies the spec- 
trum steepens, seen as a soft excess over the hard X-ray power-law 
with a featureless, approximately power-law form. The shape of 
the soft excess is similar to that observed in other objects. For ex- 
ample, Marshall et al. (2002) found the high-resolution soft X-ray 
spectrum of the ultrasoft Seyfert 1 Mrk 478 (which has a similar 
X-ray luminosity to Ton SI 80) was a featureless power-law. PKS 
0558-504 also showed a similar soft excess shape (O'Brien et al. 
2001). The Chandra HETGS observation of the variable, narrow- 
line Seyfert 1 galaxy NGC 4051 again showed the soft excess to be 
a rapidly variable continuum component. However, in this lower- 
luminosity source it showed significant spectral curvature (Collinge 
etal. 2001). 

The observed flux of Ton SI 80 showed rapid variability, 
changing by ~ 50 per cent during the observation, but showed only 
weak spectral variability. This implies that, to first order at least, 
the curved X-ray spectrum is varying as one component. There 
are indications that this may be a general result, that Seyferts with 
steep X-ray spectra show very little spectral variability (at least on 
short timescales). The ASCA monitoring of Ark 564 and Ton S 180 
showed flat RMS spectra (Turner et al. 2001c; Romano et al. 2001; 
Edelson et al. 2002), as did the recent XMM-Newton observation 
of IH 0707-495 (Boiler et al. 2002). This lack of short timescale 
spectral variability is somewhat different from the case in "normal" 
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Seyfert Is, which tend to show stronger variability at lower energies 
(e.g., Nandra et al. 1997; Markowitz & Edelson, 2001; Vaughan & 
Edelson, 2001). 

The lack of strong spectral features means that the observed 
form of the soft excess is consistent with a range of models. The 
data clearly rule out models based on blends of narrow soft X- 
ray lines (e.g.. Turner et al. 1991), but leave a variety of contin- 
uum emission mechanisms as plausible alternatives. Reprocess- 
ing of primary X-rays by an ionised accretion disc can produce 
a strong, smooth excess similar to that observed but has di fficult y 
simultaneously explaining the hard X-ray spectrum (section 3.2.3). 



The observed lack of spectral variability is also a challenge for this 
model. The highly correlated, simultaneous variability between 3- 
10 keV flux (dominated by the primary X-rays) and softer X-ray 
bands (dominated by the reprocessed emission) requires the repro- 
cessing to produce a near-perfect "reverberation" signature with a 
delay of At < 500 s (section Assuming a black hole mass of 
> 10^ Mq (from the Eddington limit) this places the reprocessor 
within <^ lOrg of the hard X-ray source. 

As first pointed out by Bechtold et al. (1987) the soft excess is 
too hot to be "bare" thermal emission from the accretion disc. The 
spectral form is consistent with emission from multiple blackbody 
components, but the derived temperatures and sizes are inconsis- 
tent with those predicted for an accretion disc (unless Ton SI 80 is 
highly super-Eddington). Allowing for Doppler and gravitational 
shifts does not significantly alter this result. 

Inverse-Compton scattering of soft photons by thermal elec- 
trons provides a more physically s atisfy ing explanation for the 
broad soft excess (sections 3.2.5 and 3.2.( ), in which case the seed 



photon source is consistent with thermal accretion disc emission. 
The harder power-law emission extending to 10 keV (and presum- 
ably beyond) can also be produced by Comptonisation, either in 
another purely thermal plasma or by non-thermal electrons in a 
plasma with a hybrid thermal/non-thermal distribution. Whether 
this is dominated by thermal or non-thermal electrons is impossi- 
ble to tell without higher energy data; fits using thermal and hybrid 
thermal/non-thermal models are comparable in the XMM-Newton 
band, where the predicted spectra both resemble power-laws. 

A problem with the models discussed above is that they have 
difficulty explaining the rapid and energy-independent variability 
of Ton SI 80 (section These constraints force the hard and soft 
X-ray producing regions to be in very close causal contact with one 
another. If at high accretion rates the disc is puffed-up due to ra- 
diation pressure, then it is possible that the separation between the 
hot corona and the optically thick accretion disc can be negligibly 
small. In this scenario, some fraction of the accretion energy is dis- 
sipated in the surface layers of the disc, producing hot electrons 
which can Comptonise the soft photons from below. The accre- 
tion disc model of Hubeny et al. (2001) demonstrates the formation 
of a hot upper layer to the accretion disc due to dissipation near 
the surface (see their figure 10), and simulations of magnetized ac- 
cretion discs by Miller & Stone (2000) suggest that even higher 
temperatures may be reached, leading to the formation of a true 
corona. Simulations of emergent disc spectra accounting for dissi- 
pation in the upper atmosphere will be presented elsewhere (Ross et 
al. in prep.), but early indications suggest that the observed strong, 
smooth soft excess can be reproduced. 

We also note the lack of a significant narrow, neutral iron emis- 
sion line at 6.4 keV. The formal (90 per cent) upper limit on the 
equivalent with of such a line is 60 eV, similar to that obtained for 
IH 0707-495 by Boiler et al. (2002). Again this contrasts with the 
situation in more normal Seyfert Is, which tend to show a narrow 



6.4 keV emission line from distant, cold material (e.g., Yaqoob, 
George & Turner, 2002). However, whether this is due to a system- 
atic difference between ultrasoft and normal Seyfert Is remains to 
be seen. 
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